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Abstract 

A problem in nonlinear water-wave propagation on the surface of an inviscid, 
stationary fluid is presented. 

The primary surface wave, suitably initiated at some radius, is taken to be a 
slowly evolving nonlinear cylindrical wave (governed by an appropriate Korteweg-de 
Vries equation); the depth is assumed to be varying in a purely radial direction. 

We consider a sech 2 profile at an initial radius (which is, following our scalings, 
rather large), and we describe the evolution as it propagates radially outwards. This 
initial profile was chosen because its evolution over constant depth is understood 
both analytically and numerically, even though it is not an exact solitary-wave 
solution of the cylindrical KdV equation. The propagation process will introduce 
reflected and re-reflected components which will also be described. The precise 
nature of these reflections is fixed by the requirements of mass conservation. 

The asymptotic results presented describe the evolution of the primary wave, 
the development of an outward shelf and also an inward (reflected) shelf. These 
results make use of specific depth variations (which were chosen to simplify the 
solution of the relevant equations), and mirror those obtained for the problem of 
1-D plane waves over variable depth, although the details here are more complex 
due to the axi-symmetry. 
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1 Introduction 



The propagation of plane solitary waves over variable depth is now well understood ([[| 
14, [H], [11], |7]). In this problem it is shown that, starting from an initial solitary- wave 



solution of the KdV equation (||12| ), the propagation over a region of varying depth 
introduces a shelf directly behind the primary wave, a left-going (reflected) 'shelf and 
another right-going 'shelf (re-reflection). All four of these components are required in 
order to ensure that the global 0(1) mass is conserved. 
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The study of the propagation of axi-symmetric (or cylindrical) waves is not, however, so 
well understood. These waves, governed by the cylindrical KdV equation 

u 1 



(as first written in |T3| , and used in their study of ion-acoustic waves), have additional 
geometric complexities which complicate their study 

To date, the work involving cylindrical solitons and the cylindrical KdV equation has been 
concerned with finding analytic solutions (0, @, H, [16], [17|, |l5j) or using numerical methods 
to investigate the evolution of given initial profiles ( [fL^H ). There have, however, been 
investigations involving additional physical features, for example, including an underlying 
shear flow ||; but nothing to the extent of those for the plane wave problem. 

The aim of this work is to draw upon the ideas and techniques used in the study of plane- 
wave propagation over variable depth, and apply them to the axi-symmetric problem. 
Some of the issues raised in || will play an important role; in fact, because of this, 
the sech 2 initial profile was used instead of the appropriate solitary-wave solution of the 
cylindrical KdV equation, in order to highlight the correspondence. 



2 Governing Equations 

We begin by introducing the equations and boundary conditions which will be used for 
this investigation. We make several assumptions about the ambient state of the fluid, 
which is taken as stretching to infinity in all horizontal directions, and lies between a free 
upper surface and a solid impermeable lower surface. We will consider an inviscid fluid, 
which implies that there will be no viscous stresses acting on the surface of the fluid, thus 
eliminating the possibility of the surface disturbances being driven by winds, for example 
(and also eliminating the existence of a boundary layer flow along the bottom). 

The fluid will be initially stationary (in its unperturbed state). In addition we assume 
the fluid to be incompressible, and we neglect the effects of surface tension. 

Written in component form, the governing equations - we use the Euler equation - for 
axi-symmetric problems are 

Du' 1 dp' Dw' I dp' / D d , d , d 

-Dt- = -pd? ] -DT = —pd^- 9 > [ Wheie Di = W +U d? +W d^ 

The continuity equation is expressed as 

1 9 , , dw' 

r' dr' dz' 

The boundary conditions are 

Dh' 

w = -jy- and p = Constant on z = h'(r',t'), 
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r'=0 



Figure 1: Configuration of the problem 

and 

db' 

w = u —— onz = o (r ). 

ar 

The equations are quoted here in dimensional form (denoted by the prime). We now 
nondimensionalise these equations by introducing the following scales: ho, a typical 
undisturbed depth; a, a typical wave amplitude; A, a typical wave length; see Fig(|l|). 
The appropriate speed scale is \fgh~o, and this is used to define a suitable time scale; we 
now obtain the familiar nondimensional equations 

u t + e (uu r + wu z ) = -p r ; 5 2 {w t + e (uw r + ww z )} = -p z , 

together with 

u 

u r H h w z = 0; 

r 

the upper boundary conditions become p = rj, w = rj t + eur] r on z = 1 + eq\ the lower 
boundary condition is w = ub r on z = b(r). In the above equations we have introduced 
p(r, z, t) to be the deviation from the undisturbed hydrostatic pressure distribution, and 
the parameters are defined by 

e = — and S = (1) 
ho A 

In addition to the above equations, we also have the equation for global mass conservation. 
If we consider the continuity equation together with the boundary conditions, and if 
undisturbed conditions exist far enough ahead of and behind the wave, then we can show 
that 

oo 

J r](r, t)rdr = constant, 
o 
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which is the conservation of mass for the surface wave in cylindrical geometry. We will 
suppose that when the wave is initiated at a given radius, it has a given profile; it therefore 
carries a known amount of mass, the above constant. 



These governing equations contain the familiar parameters for water-wave theory; we now 
proceed with the conventional transformation relevant in a far field (0) 

e 3 e 3 e 3 5 2 5 2 e 5 

This transformation accommodates the geometric decay of the surface wave as r — > oo. 
At present the topology of the bottom, b(r), changes on the same scale as the wave 
evolves, so we introduce a parameter a according to 

b(r) = B(ar), 

which will allow us to vary the scale on which the bottom changes (0). Introducing 
these into our governing equations we obtain 

U T + A(UU R + WU Z ) = -P R ; (2) 
A{W T + A(UW R + WW Z )} = -P z ; (3) 

U R + ^ + W z = 0, (4) 

with 

P = H- (5) 

W = H T + AUH R , (6) 

both on z = 1 + AH (the rescaled free surface) ,and 

W = UB R ouz = b(^r), (7) 
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where we have introduced 



e 4 



which will be the small parameter to be used in the problem, and will be, therefore, the 
scale on which the nonlinear wave evolves. 



There are then three cases of particular interest determined by the size of a in relation 
to A. These are 

(1) 5 2 a/e 2 = Act ; (2) 6 2 a/e 2 = A ; (3) 6 2 a/e 2 = a, A = 7a, 

forA^0,cr^0,7^0. Case (1) describes the situation where the depth changes on a 
scale slower than that of the evolution of the surface wave; it is this particular case which 
will be dealt with in this presentation. 
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As an aside, case (2) is where the depth changes on exactly the same scale as the wave 
evolves; case (3) is where the depth change is slow, but faster than the evolution of the 
surface wave. 

Before we proceed to investigate case (1) - Slow Depth Change, we must point out that 
in what follows, no assumption will be made regarding the relationship between the two 
parameters A and a: they will be treated as independent small parameters. 



3 Slow Depth Change 

To proceed we must introduce appropriate variables which are relevant to the various 
propagation modes which exist in the problem. The initial surface profile, the Primary 
Wave which initiates the process, propagates radially outwards; thus we need an outward 
characteristic and (slow) scale to describe its evolution; these are 

Y Y 

£ = ^ J F(y)dy - r, * = \ \ S{y)dy, Y = AaR, (8) 
i i 

where the initial wave is centred about Y = 1 i.e. R = 1/Aa (a large radius). 

The functions F(Y) and S(Y) are introduced to ensure that the characteristics and slow 
scales are consistent with the depth variation. (Note: choosing F = S = 1 recovers the 
classical characteristics and slow scales for constant depth.) 

We also anticipate that the propagation will initiate reflected components travelling in 
the opposite direction, thus we also need a suitable inward 'slow' characteristic 

Y 

F(y)dy + T, T = Aar. 



i 



Thus, any solution for the surface wave will be described (using a multiple-scale ap- 
proach) in terms of the variables (£, £,X,Y) and parameters (A, a). Introducing these 
new variables into our governing equations (|])-(|7|), we obtain 

AaU ( - U(. + A {FUU(. + ASUU X + AaFUU c + AaUU Y + WU Z } 

+FP 5 + ASP X + AaFP ( + AaP Y = 0; 
A {AaW c — Wf + A (FUWi + ASUWx + AaFUW ( + 

AaUW Y + WW z )} + P z = 0; 
AaU 
Y 



FUz + ASU X + AaFU c + AaU Y + + W Z = 0, 



with boundary conditions (on z = 1 + AH) 
P = H- 

W = AaH c - + AU {FH^ + ASH X + AaFH ( + AaH Y } 
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and 

W = -AaUD'(Y) on z=l-D(Y). 
Here we have introduced D(Y), defined as 

D(Y) = 1-B(Y), (9) 
to be the local depth, a more convenient function than B(Y). 

We now seek a solution to these equations in the form of a double asymptotic expansion 



oo oo 



Q~^^AV m g„ m , A^O, a^O, 



n=0 m=0 

where Q represents each of U, W,P,H. 

4 Results 



We proceed by constructing systems of equations at each order, A n <7 m , and imposing con- 
ditions (as appropriate) which ensure that the asymptotic expansions remain uniformly 
valid as |£| — > oo, \(\ — > oo . We present only the main results that occur at each order 
(using results from previous orders, where relevant): 



AV°: 



F(Y) = - (10) 



3 D 

2VdSH 00X + —H 00 H 00 £ + —Hoottf = 0. (11) 



3 D D 1 / 4 

2VDSH 01X + -(H 00 H 01 )s + jH 0l ^ = -2y^ {D 1 ' i Y 1 ' 2 H m ) Y . (12) 



U 1K + % + Any = 0; VDU 1K + A 11C + = 0, (13) 



AV: 



AV 2 : 



and 

2L» 1 /4yV2 (L) i/4 y i/2 /lii$)y = _ {DYHoo y) Y , (14) 
where we have introduced 

U n = Un(Z,X,Y)+U u (t,X,Y); 
H n = hn(S,X,Y)+An({,X,Y). 
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Primary Wave 




Y=l 



Figure 2: Schematic of Components in problem 



Equations ( flTD and (|12D , when suitably combined, imply 



FV£> / H^dY = Constant, 



the conservation of 'momentum' for the cylindrical KdV equation. 



In passing we observe the similarities between the above equations for the radially sym- 
metric problem and those for the plane- wave problem ; it is this similarity that suggests 
a similar approach to finding a solution. Note that, by virtue of the large-radius scaling, 
the geometric contribution in the cylindrical coordinates first appears in (0) , leaving ([LI]) 
as the classical KdV equation. The various components described by equations (|TTD (|T^) 
are represented schematically in Fig(^); see the descriptions that follow, where we give a 
brief outline of what can be deduced at each order. 



4.1 Primary Wave 



Equation (|11|) is a variable coefficient KdV equation; thus our sech 2 initial profile is an 
exact solution for general S(Y), however, in our particular problem we require S(Y) = 
d -5/2(y)y~ 2 / 3 , in order to remove nonuniformities; this, along with the requirements of 
conservation of momentum, yields a solution for the primary wave as 



= sech2 ( 2 jD 3/ 3 2 yi/3 - cX ) ) ' ( 15 ) 



where Aq is the initial amplitude (at Y = 1) of the wave. 
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4.2 Outward 'Shelf 



Now let us consider fll2f) ; we anticipate that, as this equation has distinct similarities to 
that of the plane-wave case, the solution to this equation will exhibit shelf-like properties. 

Let us introduce = £ — cX, to be a characteristic moving with the primary wave (with c 
the constant speed in (£,X) space), and let = be the peak. To demonstrate that the 
component ifoi does not approach zero near the tail of the primary wave, we construct 
an integral across a neighbourhood of = 0. This shows that 

Z 01 (T) = - ^^ (D^Y^r, 

is the amplitude of H i directly behind the primary wave. Here we have defined Y = Y(T) 
to be = 0, describing the path of the primary wave. 

From ( fl2|) we can also obtain a complete description of this wave component behind the 
primary wave. When we note that H 00 is exponentially small in its tails, we find that 
defining H 01 = D~ x l 4 Y~ x l 2 G m , then 



where 



and 



Goi ~ & ( o-0 + c j D- 2 Y- 2 ^dT 



Q(T) = -m Y 1/6 D 3/4 + A Y = Y(T),D = D(T), 



r 

T = c I D- 2 ?- 2/3 dT. 



o 



At the front (outside end) there exists a transition region (a region where transition back 
to undisturbed conditions occurs) around = 0; similarly, at the rear (inside end), a 
transition region exists around £ = 0. The detailed structure of these transitions can be 
obtained from (|T2p by considering the relevant asymptotic regions. They turn out to be, 
at the front, 

H Ql = D^Y l l & (D^Y^yL , 

where 

8_ A -5/2 /MfliA-aT^ fl _ 1 ^ _U (l Ajfl J- Mtt 



if) 



{^(taoh0-l)+ (3A 2 6 +^) sech 2 (9 
- tanh 6 sech 2 6 + § + 3A 2 # 2 ) } 



with 



(£ — cX) , s = constant; 



2j D3/2yl/3 
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-0.5 



1.46 1.47 1.48 1.49 115 1.51 



Figure 3: An example of the outward shelf behind the primary wave (Hqi against Y) 
using A = l,a = 0.01, A = 0.1 and Y = 1.5. 

and at the rear 

Hoi = D-^Y-^Q'(0) 1 1 - £ Ai{u)du} ; x = ^ . 

In Fig(||) we show a solution for the outward shelf combining the above expressions. We 
have made use of a simple depth profile, as an example {D{Y) = Y~ 2 ), and have made 
necessary choices for the parameters and constants, which allow the three asymptotic 
solutions to be suitably joined at this order of approximation. 



4.3 Mass Conservation 

As we have already mentioned, the specific nature of this problem leads to an underlying 
requirement for global conservation of mass. After the application of all of the relevant 
scalings we obtain the mass conservation condition as 

oo 

H (R, r)RdR = J (H 00 + aH 01 + AH W + AaH n + ... )RdR 

o o 

= Constant. 

Now let us consider the mass carried by the primary wave; this turns out to be 

oo 

H m RdR~D l/2 Y [ H m d(b = m n DY 2/3 , 



o 



a quantity which is clearly not constant. Thus the resulting difference from the initial 
0(1) mass, mo, must be carried by other wave components. 

We can calculate the mass carried by the Outward Shelf as 

oo 

o J H 01 (£, X, R)RdR ~ m (D 1/4 Y 1/2 - DY 2/S ) , 
o 

which satisfies the condition that the mass carried by this component is zero at Y = 1. 
We can see that the total mass carried by these two components is 



J (ifoo + crH 01 )RdR = m D 1/4 Y 1/2 + O(A), 
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Thus there must be other components carrying 0(1) mass. 



4.4 Re-reflected Shelf 



Let us now investigate the mass carried by the outward moving component of H\\. (We 
have written the non-local contributions to Hu as tin + An; see (|13),(|14).) This compo- 
nent is governed by flH|) which can be integrated to yield 

Y 

h m = -Vv4y-V2 J D-^Y-y^{DYH 00Y }dY. 

1 

For general depth profiles we are unable to obtain an explicit description for this wave 
component. However, the use of a specific form for the depth change, namely 

(/3 + 7 y2)4/3 



D(Y) 



Y 2 



(16) 



enables us to completely solve this equation and thereby elucidate some properties of this 
component. 



Calculating the mass carried by this component, Mo, we now find that 

Y 

M = m 



2 7 (/3 + 7F 2 ) 2 / 3 4(/3 + 7F 2 ) 4 / 3 



y4/3 
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+ (--2 T )F(^ + 7 y 2 ^ 1 ; 



dY. 



4.5 Reflected Shelf 

This component is governed by the pair of linear homogeneous equations ([13]). In order 
that we have conserved 0(1) mass, this component must be carrying the mass 

m - m D 1/4 Y 1/2 - M , 

and this requirement uniquely determines An. The details are involved, but readily 
accessible (at least for suitable D(Y)). 

The full solutions for both the reflected and the re-reflected shelves, for our specific choice 
for the depth profile flT6|), are not presented here; but they can be seen, along with the 
description of the transitions regions for the re- reflected shelf, in || . 



5 Conclusions 



In this work we have presented a description of the wave components in cylindrical geom- 
etry that carry 0(1) mass. It has, however, not been possible to present all the details for 
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general depth profiles (which would have been preferable). Nevertheless, we have shown 
the development of the method of solution, and from this can obtain some particular 
results. 

Thus, for example, we find that if the depth profile takes certain (very simple) profiles, 
e.g. D(Y) oc F~ 2 / 9 , then one of the wave components, the new outward (re-reflected) 
shelf, will be absent. 

One fundamental question which still has to be answered, concerns the behaviour of the 
inward (reflected) wave component. It is obvious, by the very nature of the problem, 
that any component propagating inwards will eventually reach the centre. There will 
undoubtably be some form of singularity at the centre, but there may be some choice for 
the depth profile in the neighbourhood of r = which will result in a particularly simple 
structure. The details of this investigation are still being developed, and may provide an 
impetus for further study. 
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